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I We consider the formation of chiral density waves in Quarkyonic matter, which is a 

phase where cold, dense quarks experience confining forces. We model confinement 
following Gribov and Zwanziger, taking the gluon propagator, in Coulomb gauge 
and momentum space, as ~ We assume that the number of colors, iVc, is 

large, and that the quark chemical potential, fi, is much larger than renormalization 
mass scale, Aqcd- To leading order in 1/Nc and AqcDZ/ti, a gauge theory with A'f 
Q flavors of massless quarks in 3 + 1 dimensions naturally reduces to a gauge theory 

00 in 1 + 1 dimensions, with an enlarged flavor symmetry of SU{2Ni). Through an 

^ anomalous chiral rotation, in two dimensions a Fermi sea of massless quarks maps 

directly onto the corresponding theory in vacuum. A chiral condensate forms locally, 
and varies with the spatial position, z, as {'ip exp(2ifj,z'y^^^)ip) . Following Schon and 
O Thies, we term this two dimensional pion condensate a (Quarkyonic) chiral spiral. 

>■ Massive quarks also exhibit chiral spirals, with the magnitude of the oscillations 

decreasing smoothly with increasing mass. The power law correlations of the Wess- 
;J] Zumino-Novikov-Witten model in 1 + 1 dimensions then generate strong infrared 

^ effects in 3 + 1 dimensions. 



1 Introduction 



The phases of Quantum Chromodynamics (QCD) at nonzero temperature and 
density are a subject of continuing interest. While numerical simulations on 
the lattice can be of use at nonzero temperature when the quark density is 
small, standard Monte Carlo techniques are not of use in cold, dense quark 
matter. 
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One expansion which is of utihty is to expand in the hmit of a large number of 
colors [Tll2]|Sllill3] ■ For cold, dense quark matter — quarks in the fundamental 
representation, coupled to an SU (Nc) gauge theory — this gives a "Quarky- 
onic" phase [6.7.8.9J . Keeping the quark chemical potential, /x, of order one as 
the number of color Nc oo, the free energy for this phase is dominated by 
that of quarks. Nonetheless, excitations near the Fermi surface are confined, 
perhaps baryonic, whence the name. 

In this paper, we consider chiral symmetry breaking in Quarkyonic matter. 
We consider a phenomenological model for confinement, taking the timelike 
component of the gluon propagator to be D^'^ ~ This is valid in 

Coulomb gauge, for a spatial momentum p, and corresponds to a potential 
which rises linearly in coordinate space. Such a propagator was originally 
suggested by Gribov [TD] and Zwanziger [ll]. To use such a propagator in 
cold, dense, quark matter, it is necessary to assume that gluons are insensitive 
to screening by quarks. For this to be true, the number of flavors, A^f , must be 
<C Nc, and the chemical potential must satisfy fi <^ N^^'^Aqqy), where Aqcd 
is the renormalization mass scale of QCD [6j|. 

Chiral symmetry breaking in such a model has been studied by Glozman and 
Wagenbrunn [8] and by Guo and Szczepaniak [9], for values of ~ Aqcd- 
We work in the extreme Quarkyonic limit, /i ^ Aqcd, so that the effects of 
chiral symmetry breaking in vacuum can be ignored. It is possible for chiral 
symmetry breaking to occur at large /i, since we are, by assumption, in a con- 
fined regime. Of course there is no guarantee that our results apply to QCD, 
where Nc = Nf = 3; nevertheless, there is certainly some range of Nc, N^, and 
/I, where it does. If applicable to QCD, our results are of interest to inter- 
mediate densities, where both conventional nuclear physics and perturbative 
treatments fail. Notably, this may include the astrophysics of neutron stars. 

In vacuum, chiral symmetry breaking occurs through the pairing of a left 
handed quark with a right handed anti-quark, {ip^ipL) 7^ 0, and vice versa. 
This condensate is, of course, spatially uniform, so that the spontaneous break- 
ing of chiral symmetry does not disturb the Lorentz invariance of the vacuum. 

Now consider the effects of a Fermi sea, where there is a net excess of quarks 
over anti-quarks. The analogy of the usual condensate is illustrated in Fig. 
[l] Energetically, it costs essentially zero energy to excite a quark right at the 
edge of the Fermi sea. On the other hand, it costs at least ~ 2/i to pull an 
anti-quark out from deep in the Dirac sea. (Remember that we assume that 
/i is very large.) Thus the usual condensate can not be formed spontaneously, 
and anti-quarks will not enter into our analysis henceforth. 

There are numerous features which are not captured by the illustration in 
Fig. [TJ We really should draw not one, but two Fermi seas: one for left handed 
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Fig. 1. Quark anti-quark condensate in the presence of a Fermi sea. The energy of 
the pair is E'tot — 2/i, while its momentum is Ptot — 0. 

quarks, and one for right handed quarks. To avoid unnecessary dupHcation, 
instead we assume that the quark, denoted by a filled circle, is always left 
handed, and that the anti-quark, denoted by an open circle, is right handed. 
The quark and anti-quark are also assumed to have the same color, so that 
any condensate is ~ A^c, and survives in the limit of large Nc. If the quark 
has momentum p, then the anti-quark, formed by removing a quark with 
momentum p from the Dirac sea, has momentum —p. Thus the quark anti- 
quark pair has no net momentum, and this condensate is spatially uniform, 
as in vacuum. 



In the presence of a Fermi sea, though, it is also possible for chiral symmetry 
to be broken by pairing, say, a (left-handed) quark and a (right-handed) quark 
hole. If both the quark and the quark hole are near the edge of the Fermi sur- 
face, then it costs little energy to excite them, and the energetic penalty paid 
to excite an antiquark can be avoided. The natural analogy to the condensate 
in vacuum is illustrated in Fig. [2} pairing a quark with momentum p, and a 
hole, formed by removing a quark with momentum p from the Fermi surface. 
The momentum of the hole is then —p, so the quark-hole pair has no net mo- 
mentum, and is spatially constant. In condensed matter physics, an excitation 
as in Fig. |2]is known as an exciton. Naively, we might expect that excitons are 
suppressed, since the relative momentum between the particle and the hole, 
2p, is large. 

However, this is not the only way for quarks and their holes to break the chiral 
symmetry. Consider pairing a (left handed) quark, with momentum p, and the 
hole formed by removing a (right handed) quark with the opposite momentum, 
—p, from the Fermi sea. The quark hole then has the same momentum as the 
quark, +p, so that the resulting condensate is not uniform, and has a net 
momentum +2p; this is, it varies as ~ exp{2i^z), where z is the direction 
along which the pair moves, p = pz. Such condensates do not occur in vacuum, 
where they would imply the spontaneous breaking of rotational symmetry. In 
condensed matter physics, though, such non-uniform condensates are common, 
and known as density waves [H]; this is then a chiral density wave. Note that 
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Fig. 2. Exciton pairing between a quark and a quark hole; E'tot — Ptot — 0. 

the relative momentum between the quark and its hole is small, so such a 
condensate may be favored. 




Fig. 3. Pairing between a quark and a quark hole, which generates a chiral density 
wave; E^tot - 0, Ptot - 

In this paper we show that in the Gribov-Zwanziger model, that the exciton 
pairing of Fig. [2] is not generated, but that the chiral density wave of Fig. |3]is. 
Again, this is familiar from systems in condensed matter: typically excitons are 
only created dynamically as resonances, such as by the absorption of light, and 
usually do not condense. Density waves are common, especially for systems in 
1 + 1 dimensions [Tl|. We will investigate all Dirac and flavor structures, and 
show which types are preferred. 

For completeness, we illustrate the pairing between two quarks which leads to 
color superconductivity in Fig. |4} This is pairing between a quark at one edge 
of the Fermi surface, with momentum +p, and another quark at the other edge, 
with momentum —p. Since pairing is between two quarks, the condensate has 
no net momentum and is spatially uniform. For this reason, pairing can occur 
over the entire Fermi surface, in a spatially symmetric state. 

So far, we have not emphasized the Nc and iVf dependence of pairing, which is 
not captured by the illustrations in Figs. [l]-|4} The pairing in Figs. [2] and |3] is 
between a quark and a quark hole of the same color, so the condensate is ~ N^. 
Further, to the extent that A^f -C N^., the condensate is rather insensitive to 
Af. In contrast, the diquark pairing of color superconductivity depends upon 
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Nc and A^f in an essential way. Fermi statistics greatly constrains the pairing 
between two quarks (or two quark holes): it is always anti-symmetric in color, 
so there are strong relations between the spatial wavefunction, flavor, and 
chirality. For instance, in case of A^f = 2 and A^c = 3, spatially symmetric 
condensates form by anti-symmetrizing in flavor; this condensate pairs quarks 
of the same chirality together, and so does not break the chiral symmetry. On 
the other hand, for A^f = 3 and A"c = 3, the preferred condensate does break 
the chiral symmetry, through color- flavor locking [12] . 

For more than three colors, the gaps for color superconductivity depend sen- 
sitively upon which representation one assumes the quarks to lie in. If the 
quarks are in the fundamental representation, then since the pairing for color 
superconductivity is anti-symmetric in the colors of the two quarks, the gap is 
not a color singlet, and is suppressed at large A'c. It is also possible, however, 
to generalize QCD by letting the quarks lie in the two-index, anti-symmetric 
representation of color [13j. This limit is rather different from that which we 
consider in this paper. There are ~ A"^ quarks in this limit, so that gluons are 
affected the quarks, and there is no Quarkyonic phase. This is like taking the 
number of flavors, A^f , to grow with N^. In such a limit color superconductivity 
is not suppressed at large A'c. 

It is not clear which of these two limits is most like QCD, with three colors 
and three light flavors. We suggest that it is useful to consider all possible 
hmits, and to see what qualitative conclusions might be tested in QCD. 
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Fig. 4. Pairing between two quarks which generates color superconductivity; 

^tot ^ Ptot ^ 0. 

If a channel for color superconductivity exists, then Cooper pairs will form for 
arbitrarily weak coupling. Thus color superconductivity is always the domi- 
nant pairing mechanism at asymptotically large chemical potential. The es- 
sential question is then, how large does the chemical potential have to be for 
color superconductivity to win out over other pairing mechanisms, such as 
chiral density waves? 

That chiral density waves [14] dominate at large A"c was first demonstrated 
by Deryagin, Grigoriev, and Rubakov [TS|16f 17II18|19|20II21|22|23I|24] . Using a 
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perturbative gluon propagator, ~ in momentum space, they find that 
chiral density waves form, with a condensate 



in magnitude, where as = g'^/'iir is the QCD fine structure constant, measured 
at a scale /i. Implicitly, the computation assumes that dense quarks form a 
Fermi liquid, so that pairing is from quarks (and holes) within Apert of the 
edge of the Fermi sea. 

In contrast, Quarkyonic matter is not a Fermi liquid because of confinement 
of quarks. Low energy excitations, within ~ Aqcd of the edge of the Fermi 
sea, interact not through the perturbative gluon propagator, but through the 
Gribov-Zwanziger form, ~ ^)^. A chiral density wave forms, with the same 
Dirac structure as in perturbation theory. Because the relevant momenta for 
pairing is controlled by a confining potential, the Quarkyonic condensate for 
chiral density waves is inevitably ~ Aqcd in magnitude. 

Our analysis applies in a regime of intermediate /z, where the Quarkyonic gap is 
greater in magnitude than the perturbative gap. we find that for Apert < Aqcd 
whenever > 0.12, using Aqcd ~ fiexp{—6n/{llNcas)). For the coupling 
in QCD, this means that < 80 GeV. At larger values of /i, Apert > Aqcd- 
As we discuss at the end of Sec. 2.1, this is a difficult regime to treat, as the 
effects of both perturbative interactions, and confinement, must be included. 

It is known that for a Fermi liquid of dense quarks, chiral density waves 
lose out to color superconductivity except for very large values of the color, 
Nc > 1000 iVf P^IITT] . This is because the chiral density waves generated by 
perturbative interactions are very sensitive to screening by dynamical quarks. 

In sharp contrast, we expect that Quarkyonic chiral density waves are much 
less sensitive to screening by dynamical quarks. At large Nc, dynamical quarks 
do not affect a Quarkyonic phase until asymptotically large values of the 
chemical potential, /i ~ N^^'^Aqqy). This power of N^^'^ follows either from 
considering the Debye mass, as in Ref. [B], or the free energy, as discussed in 
Appendix |Aj As long as there is a Quarkyonic phase, and Apert < Aqcd, we 
expect that Quarkyonic chiral density waves dominate. 

The outline of the paper is as follows. For quark and quark hole excitations 
near the Fermi surface, where the magnitude of the transverse momentum 
\p± I ^ the theory in 3 + 1 dimensions reduces to an effective model in 1 + 1 
dimensions [TipB] . In Sec. [2] we show that in the Gribov-Zwanziger model. 




(1) 
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this implies that in the gluon propagator, we can integrate over f)± to obtain 

This is the gluon propagator in 1 + 1 dimensions, so our effective theory is just 
QCD in 1 + 1 dimensions. We discuss how in the Gribov-Zwanziger model, it 
is necessary for /x ^ Aqcd for this reduction to hold. 

In Sec. [3] we consider how quantum numbers in 3 + 1 dimensions map onto 
those in 1 + 1 dimensions. Starting with left and right handed massless quarks 
in 3 + 1 dimensions, we find that the reduced model has a doubled flavor 
symmetry: A^f flavors in 3 + 1 dimensions becomes an SU{2Nf) symmetry 
in 1 + 1 dimensions. This extended symmetry follows immediately from the 
analysis of Shuster and Son (TU], and is very much like the doubling of flavor 
symmetry which occurs for heavy quarks. 

In Sec. |4]we show how in 1 + 1 dimensions, through an anomalous redefinition 
of the quark fields, a theory at nonzero chemical potential can be mapped 
onto the corresponding theory in vacuum. The net quark number, present in 
the theory at /i 7^ 0, is generated by the axial anomaly, as shown in Appendix 
[Bj The mapping can then be reversed: knowing results for a gauge theory in 
1 + 1 dimensions in vacuum, one can read off what happens in a Fermi sea, 
/i 7^ 0. We show that a constant condensate in vacuum, (ip'ip) 7^ 0, produces a 
condensate {ip exp{2inzr^)ip) 7^ 0, where z is the spatial coordinate, and the 
Dirac matrix in 1 + 1 dimensions. In the reduced model, the two dimensional 
= 7o7z, where 70 and 72 are Dirac matrices in 3 + 1 dimensions. 

This type of spatially dependent condensate at /i 7^ is familiar from soluble 
models in 1 + 1 dimensions [18], where Schon and Thies termed it a "chiral 
spiral" [19j. A chiral spiral was also found by Bringoltz, in his numerical 
analysis of heavy quarks at nonzero density in QCD in 1 + 1 dimensions [20] . 
We thus term our solution a Quarkyonic Chiral Spiral (QCS). 

While we concentrate on massless quarks, in Sec. [5] we discuss how massive 
quarks also exhibit chiral spirals. This is because even for massive quarks, 
excitations about the Fermi surface are gapless at tree level. 

Effective theories about a QCS are considered in Sec. |6} Using non-Abelian 
bosonization, the reduced model in 1+1 dimensions reduces to a Wess-Zumino- 
Novikov-Witten (WZNW) model [25|26f27j . This model has long range corre- 
lations, which should also produce long range correlations in 3 + 1 dimensions. 

We conclude in Sec. [7] about whether QCS's may be be relevant for cold, 
dense QCD. We note that QCS's are closely related to pion condensation [22] . 
They are not identical, because a pion condensate is a chiral spiral in the four 
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dimensional 75, while for a QCS, it is in the two dimensional F^. We also note 
that if pionic chiral spirals occur, then probably so do kaonic chiral spirals. As 
a spatially varying condensate, a kaonic chiral spiral differs for the spatially 
constant condensate of kaon condensation [23] . Chiral density waves also arise 
in the Sakai-Sugimoto model [2T] . 

The physics of QCS's should be especially rich, however, since it includes 
the spontaneous breaking of translational and rotational symmetries, and a 
plethora of light modes. Such phenomenon should have direct implications for 
observations of neutron/Quarkyonic stars. 



2 Dimensional Reduction: Self-Consistent Equations 



2. 1 Reduction of the Schwinger- Dyson Equation 



We start by considering the Schwinger- Dyson equation for the quark self- 
energy. In a Quarkyonic phase at large A^'c, the gluons are unaffected by the 
quarks, so that corrections to the gluon self-energy, and vertices, can be ne- 
glected. For the quark self-energy we take the sum of rainbow diagrams, 

m = - / K'^ip - ^) (7.^^) s{k) ii.ts) , (3) 

where is the quark self-energy. At a nonzero chemical potential /i, the 
dressed quark propagator, S{k), is 

where S is determined self-consistently through the integral equation, Eq. (|3|. 
We work in Euchdean spacetime, k = {ki,k). In the limit that ^ Aqcd, 
we neglect chiral symmetry effects as in vacuum, as illustrated in Fig. [Tj such 
effects have been considered at /i 7^ by Refs. [SI'S]. Consequently, we neglect 
terms ~ 1 in the quark self-energy and propagator. 

For the gluon propagator, Dn^{k), we take the Gribov-Zwanziger form, 

Dff(k) = -^x^S^'' ; D'^ = D^^ = 0. (5) 
Cf (A;2)2 

At the outset, we stress that we are dealing with a model of confinement. While 
the gluon propagator, and vertices, are unaffected by quark loops, there is no 
fundamental justification in taking the vertices to be the same as the bare 
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ones, nor in taking the gluon propagator of Eq. The gluon propagator 
involves a parameter, a, which is the string tension, and has dimensions of 
mass squared, a ~ Aq^j-,. Numerical factors Svr and Cp = {N^ — l)/2Nc are 
multiplied to reproduce a correct linear potential for the color singlet channel. 

We stress that the propagator in Eq. ([s]) is vahd only for small momenta, for 
k < Aqcd- For larger momenta, one should use the usual gluon propagator of 
perturbation theory, ~ For excitations near the edge of the Fermi sea, 

though, we can neglect the perturbative part of the propagator. This differs, 
for example, from the computation of the free energy [B]. That is dominated by 
momenta transfers within the entire Fermi sea, ~ /i, for which the perturbative 
gluon (and quark) propagators should be used. There are contributions to the 
free energy from momenta ~ Aqcd, but this are small, powers of ~ (Aqcd//^)^ 
times the perturbative terms [6l 



After summing over the color indices, Eq. ([s]) becomes 

^(p) = [ -^I^ ^4 S{k) 74 ; (6) 

normalizing the generators as ti^f^t^) = 6^^ /2. Since the right hand side of 
this equation is independent of p^, so is the quark self-energy, ^(p). 

The Schwinger-Dyson equation simplifies considerably if we consider only exci- 
tations near the edge of a Fermi sea. For a free massless quark with momentum 
p'^ = {p4:,Pz,Py,Px) = {p4.,Pz,P±), in a Fermi sea its mass shell is given by 



m = \/pI+p±- l^- (7) 

Assuming that quark is along the z direction, so that \pz\ = fi + Spz, 

-'2 

ipi^5p, + ^ + ... . (8) 

Thus, as is well known JT^, in a Fermi sea the dispersion relation linearizes 
in 6pz, allowing us to neglect the effects of the transverse momenta for the 
quarks. 

Since this is the mass shell for a free quark, it neglects the effects of the quark 
self-energy, S. We expect, though, that for quarks and quark holes near the 
Fermi surface, including the quark self-energy does not affect the suppression 
of fluctuations in p±. For the Gribov-Zwanziger potential, the natural scale 
for the transverse momenta is ~ Aqcd- Thus we can neglect the quark 
transverse momenta in the extreme Quarkyonic limit, where /i ^ Aqcd- What 
happens when /i ~ Aqcd is a difficult problem which we do not address here. 

In the Schwinger-Dyson equation, the dominant contribution from such an 
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infrared singular gluon propagator is when the gluon momentum is smalL 
This constrains the internal and external momenta of the quarks to be near 
one another, k ^ p. Neglecting the transverse momenta of the quark, we 
need only consider the two components of the momenta along the light cone, 
kz and k^. The Schwinger-Dyson equation thus reduces to 

m ^ / 74 S{k,, kz, Ox) 74— — ^ . (9) 

Since the transverse momentum k± enters only through the gluon propagator, 
we can now integrate it out, 

[ ^ ^ Nc9Ib 1 . ^^ 2 

(10) 

where we define a two dimensional gauge coupling constant (720, and neglect 
The reduced Schwinger-Dyson equation then becomes 

Up^,p., Ox) ^ ^ / 74 s{k„ kz, Ox) 74 Tvr^ • (11) 

2 J (27r)^ (kz-pz) 



This integral equation corresponds to QCD in 1 + 1 dimensions. In axial gauge, 
Az = 0, the two dimensional action for gluons reduces to a free term, 

^tr Gl^ = tr {dzA,f . (12) 

In the limit of large N^, the model in 1 + 1 dimensions obviously gives the 



integral equation of Eq. (11). 



Of course it is necessary to take some care in the reduction of the Dirac 
matrices from 3 + 1 to 1 + 1 dimensions. We address this in Sec. [Sj 

A more realistic model for the gluon propagator than Eq. ([s]) is to take a sum 
of the Gribov-Zwanziger term, ~ l/(fc^)^, plus a perturbative piece, ~ 1/fc^. 
The confining term is valid for momenta < Aqcd! the perturbative term, for 
momenta > Aqcd- The analysis goes through as above. The dependence of the 
quark propagators on the transverse momenta can be neglected, so one is left 
with an integral of the gluon propagator with respect to k±. Integration over 
the Gribov-Zwanziger propagator gives ~ /{kz —PzY, while the integral over 
the perturbative piece generates a logarithm P^|T7| . ~ CpQ^ log(/^^/ (kz—pzY)- 

With only a perturbative propagator, cr = 0, the analysis is as follows p!^16fl7j : 
Assuming that the is a constant gap, Apert, the dimensional reduction 
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can apply to the region of quark momenta, A^^^^f fi ^ \Skz\ ^ Apcrt (Skz is 
measured from the Fermi momentum). Integrating over of a quark prop- 
agator yields a factor of the inverse energy, l/lSk"^ + Apgj.^|^/^. The form 
of self-consistent equation has a similar structure as the BCS gap equation 
for the constant gap. A dominant contribution to the integral comes from 
the inverse energy part, which is sensitive to the gap (an absence of gap 
yields infrared divergence in the integrand. This is nothing but the instabil- 
ity of the Fermi surface.) Integration over k^ at soft momentum region gives 
a logarithm term, ~ log(/i/Apert)- Combining it with the logarithm of the 
gluon propagator, one find the squared logarithmic form of the gap equation, 
Apcrt ~ fi'^A'c log^(;u/Apcrt) Apcrt, with the self-consistent solution of Eq. (!]). 

When we simply add the confining interaction to the self-consistent equation, 
its nature strongly changes. The self-energy can be no longer identified as 
the gap from the symmetry breaking, because the single-excitation-energy 
gap due to the confinement potential is large or divergent. The integral over 
kz in the gap equation is / dk^l/lSk'^ + S^l^''^ x a/{5kz — Sp^Y ~ cr/ J^m. x 
f{Spz), where f{Spz) is a some regular function and Air is an infrared cutoff. 
The infrared cutoff Air for the gap equation for the single particle spectrum 
should be identified with the Debye mass, so that the gap itself is of order 
^/NcAQQJ) / {^/g^^^fi) . To see this imagine we try to ionize a hadronic state into 
its constituent single particle quark excitations. This can only happen when 
the separation between the constituents is of order the Debye screening length, 
and costs an energy cri?Debye = V^-^qcd / iV d'^ ^cfJ') ■ This provides a gauge 
invariant definition for the single particle excitation gap. It also shows that the 
gap is large compared to the confinement scale in the Quarkyonic Phase. There 
is no weakly coupled solution for the gap equation when in this phase since, if 
one goes back to the derivation of the perturbative contribution, the derivation 
of the logarithmic terms is no longer valid. The dominant contribution for the 
perturbative piece in arises for momenta greater than the gap itself. 

It is important to understand that the analysis here concerns only the single 
particle excitation spectrum. The issue of chiral condensation is not directly 
related. We have provided a self-consistent derivation of the chiral condensate 
in other parts of this paper, and the condensate arises from non-perturbative 
effects. In a self-consistent analysis of course one cannot rule out the possibility 
that there is another solution arising from a different kinematic region. 

As discussed in the Introduction, we consider only a region of intermediate 
/i, where the Apcrt < Aqcd- In this region, it is safe to include only the 
confining propagator, and neglect the perturbative term. When Apcrt > Aqcd, 
we enter a more complicated regime, as the large difference in momentum 
scales required by the perturbative dimensional reduction, between A'^^^^/fi 
and Apert, is manifestly affected by confinement. Confinement can probably 
be neglected when Aqcd < ^pert/ f^y this only occurs when ~ 0.03, 
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which is an astronomically large value of fi, ~ lO^^GeV. 



2.2 Reduction of the Bethe-Salpeter Equation 



In this subsection we show how the reduction of the Schwinger- Dyson equation 
for the quark propagator applies as well to the Bethe-Salpeter equation for 
mesonic wave functions. Henceforth we work in Minkowski spacetime, k = 
{kQ,k). We consider the homogeneous equation, which in large is of the 
ladder type, 



ab 



d^k 

(2^ 



S{k-P){rtA) ^{P;k) {YtB)S{k+P) 



(13) 

a and b denote color indices, a and /3 spinor indices. \l/(P;g) is a bound 
state wave function, where 2P is the total momentum, q the relative momen- 
tum. The dressed quark propagator, S{q), is the self-consistent solution to the 
Schwinger-Dyson equation, Eq. ([TT]). The Bethe-Salpeter equation includes 
both singlet and adjoint channels, = S""^"^ + t'^'$A', for the singlet channel, 
this reduces to {Dqq = —1)44) 



*(P;g), 



d^k r 



(2vr) 



S{k-P) 7°^(P;A;)7° S{k + P) 



°/3((A;-g) 



2^2 



• (14) 



As in the previous subsection, we concentrate on very soft gluons, k ^ q. 
We assume that the both quarks, with momenta k — P and k + P are close 
to the Fermi momentum, \pf\ = /i. We assume that the quark and quark 
hole pair have a total momentum along the z-direction, 2P = 2Pz, and that 
the transverse momentum of the wavefunction can be neglected. We consider 
pairing energies that are both of the exciton type, P^ ^ and Ifc^l ~ pp, and 
for a chiral density wave, |Pj| ~ p_F and k^ — 0. 

Neglecting the transverse momenta of the quark and quark hole, as we did for 
the quark mass shell in Eq. ([s]), we then introduce a wavefuntion, Lf, which is 
obtained by averaging \1/(P; q) over the relative, transverse momenta, q±: 



(fq±_ 
f27r)2 



^(P;go,g^,gl) 



(15) 



The Bethe-Salpeter equation satisfied by this averaged wavefunction is 



i2n 



- S{q - P)a^S{q + P)si3 



d^k 



7°^(P; A;)7° 



Svra 
75 ((A;-g)2)2' 



(16) 
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We take the momenta of the quarks to be near the Fermi surface, q±P 
(pi? ± 6pz, ±6p±). Thus we introduce two dimensional quark propagators, 



S{P±q) = S{Po±qo,Pz±qz,0_ 
The Bethe-Salpeter equation becomes 



(17) 



z)al3 S{q- P)a^S{q + P)sf} 



X 



(27r) 



7°^(P; k)i 



d^qi. 



Ma 



-rs J (27r)2 ((fc_g')2)2 • 



The wave function on the right hand side depends only upon P and k, but 
not upon q, and so we can integrate k) with respect to k±. That implies 
that as for the quark self-energy, we can integrate over q±, leaving 



^{P; (io,q. 



2 

dkodk 



Siq-P)^^Siq + P) 
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7V(P; fc)7° 



(19) 



T"' {kz - qzY 



This is the same form as the Bethe-Salpeter equation for QCD in two dimen- 
sions [1], in = gauge. 



3 The Effective Lagrangian in 1 + 1 Dimensions 



In the previous Section we demonstrated how dimensional reduction, from 
3 -|- 1 to 1 -f- 1 dimensions, occurs for both the quark self-energy and the Bethe- 
Salpeter wave function. In this section we show how this arises at the level 
of effective Lagrangians. Our discussion elaborates upon that by Shuster and 
Son [16J. 

Start with the free quark Lagrangian in 3-1-1 dimensions. Since we concentrate 
on quarks near the Fermi surface, we assume that the spatial momentum is 
along the z direction, 

<f™'^^ = ^(V^7%V^ + V^7^5.^)- (20) 
For the Dirac matrices, we take 



" 


1 ■ 




" -(j^ ' 




" 1 " 


1 





, 7^ = 





, 7^ = 


-1 



(21) 



are the Pauli matrices, and 1 the unit matrix in two dimensions. We also 
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define spin matrices as 



a' 
a' 



(22) 



e*-^'' is the totally antisymmetric tensor, e^^^ = 1. If transverse momenta can 
be neglected, the reduced Lagrangian in 1 + 1 dimensions has an extended 
symmetry, which is related to the spin that quarks carry in 3 + 1 dimensions. 



We introduce projectors for the quark fields 

1 ± 



/ 1^75 , , 



L± 



1 ±7^7 



(23) 



i/jr and ipL are right and left handed fields, eigenstates of chirality. We also in- 
troduce projectors for spin [28] along the z direction, (l±E^)/2; for eigenstates 
of chirality, this equals (1 ± 7^7^ )/2. 



The usual chiral basis is to take 



R+ 



(24) 



but for this problem, this is rather inconvenient. To see this, we write the free 
Lagrangian in terms of 4'r± and 4'l±, 



',t /a , A , „/,t 



(25) 



The fields and ipi- have spin up along the +z direction, while ipL^ and 
ipR- have spin along the —z direction. In terms of two dimensions, for positive 
energy ip^^ and ipL+ are right moving fields, while ipi- and ipR- are left moving 
fields. 

Because there is no spin in two dimensions, we only need two component 
spinors. Thus from one four component spinor in 3 + 1 dimensions we obtain 
two types of two component spinors in 1 + 1 dimensions. 



i^R+ 



IpR- 



(26) 



The spinors (p^ and are eigenstate of spin along the z direction, and act 
as two "flavors" in 1 + 1 dimensions. This is valid to the extent that we can 
neglect the transverse momenta, which couple to the flavor breaking matrix, 
7±- 



We can then combine ip-^ and ip_i_ into one four component spinor, $, that is 
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related to the original quark field, by a unitary transform, U : 



.n. 



u 



1 

















1 














1 





1 









(27) 



For the Dirac matrices in two dimensions, V^, jji — 0, we make the obvious 
choice. 



In two dimensions. 



-la 



(28) 



(29) 



While this identity is trivial mathematically, it plays an important role in the 
next section. The conjugate is defined as $ = (<^-|-, (p^) = {ipp V^|)r°. 

It is straightforward to rewrite quark bilinears in 3 + 1 dimensions in terms of 
spinors in 1 + 1 dimensions. The operators appearing in the action. 



$t $ ^ $ r° $ ; 

$t r5 ^ = $ $ , (30) 



directly map from four to two dimensions. Because of these identities, we see 
that in the presence of gauge fields, and a nonzero chemical potential, that a 
gauge theory in 3 + 1 dimensions maps into one in 1 + 1 dimensions. 



/'2d 



ir^ (9^ + i^2dA^) + /^r° 



(31) 



Fermions in two dimensions only require two components. Thus a single, four 
component spinor in four dimensions becomes two types of two components 
spinors in two dimensions: ■0, or equivalently $, become (fi, where i — 'I, I 
is the flavor index in two dimensions, generated dynamically by dimensional 
reduction. 



We note that the mass term reduces similarly, 

= A+^L- + ^L^R+ + (i? O L) = $ $ . 



(32) 



One can write the complete dictionary to go from quark bilinears in 3 + 1 
dimensions to those in 1 + 1 dimensions. We introduce matrices for the two 



r/ = (rl,r2,r3) = C/(7°E^7°S^E3)C/t 





and 








' 


1 




" -il 




" 1 " 


1 





5 


il 




-1_ 



(33) 
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For operators which are diagonal in flavor, 





r'' 



$ = 7/; 7° f/t 



where = 1, r°, or 



ror"^ 



(34) 



For operators which are not diagonal in flavor, we note that the flavors matrices 
Tf do not mix right and left moving components, so and Tf commute with 
each other. Thus operators which are not flavor singlets transform as 



F" 
F" 



$ = 7/; 7° f/t 



pOpa Q 

F^F*^ 



(35) 



The complete list of mapping for quark bilinears is given in Table [T] 
Table 1 

Transformation between quark bilinears in 1 + 1 and 3 + 1 dimensions. 





1 


r5 






1 


1 




7° 




Tl 




-i7^ 




— i7°7^ 


T2 


-7^7^ 


27^ 


^'^^^"^'^ 




7"3 




7^ 







We can also use these results to construct the relevant effective Lagrangian for 
a model in which the gluon propagator is a sum of a Gribov-Zwanziger term, 
Eq. ([5]), and a perturbative term. Integration over the former gives QCD in two 
dimensions, while integration over the perturbative gluon propagator give a 
non-Abelian Thirring model [16]. Thus the general effective model is a gauged, 
non-Abelian Thirring model. As noted in the Introduction, we consider only 
a region of intermediate /i, where the effects of the Thirring model can be 
neglected. 



4 Mapping a Fermi Sea of Massless Quarks onto the Vacuum 



In two spacetime dimensions, for massless quarks one can directly map the 
theory at /i 7^ onto that in vacuum, /i = 0. This has been noted before, 
in various guises, in the hterature before [18], especially by Christiansen and 
Schaposnik [18]. Hopefully our discussion adds clarity. 

Consider the following transformation of the quark fields: 

$ = exp (-i /i z F^) $' . (36) 
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This transformation is defined to be the same for each of the two dimen- 
sional "flavors", if^ and 1^4,, and so preserves the flavor symmetry. Under this 
transformation, the Lagrangian becomes 



That is, by redefining the quark fields, we have completely eliminated the 
chemical potential; one has transformed the theory from one in the presence 
of a Fermi sea to that in vacuum. This happens because when T^dz acts upon 



exp(z/i2;r^), it equals /iF^F^, which by Eq. (29), equals — /iF°, and so cancels 



the term for the quark chemical potential in the original Lagrangian. 

As we discuss in Sec. |5| this is special to massless quarks, and does not hold for 
massive quarks. It also does not hold in higher dimensions, where the effects 
of transverse fluctuations obviate any such correspondence. 

We then have a quandry: there is a nonzero density of quarks in the original 
theory, with 7^ 0. The vacuum has no such density, so where did it go? 



The answer is that the transformation of Eq. (36) is anomalous, involving the 
Dirac matrix F^. 

One can show that the correct quark density is given, precisely, by the anomaly. 
There are many ways of doing the calculation; in Appendix |B] we give the 
computation based upon operator regularization. The computation also shows 
that the only quark bilinear to receive an anomalous contribution is that for 
quark number. For other operators, the transformation from $ to $' can be 
computed naively. 

The most interesting transformation is for the chiral condensate. Using Eq. 



(32), we write the a chiral condensate for $, in terms of 

= cos(2;Uz) - zsin(2/i2) $F^$ . (38) 

Assume that there is chiral symmetry breaking in vacuum, so that ($'$') 7^ 0. 
Actually, in 1 + 1 dimensions fluctuations disorder the system, and only leave 
quasi long range order [H]. Since this is due to fluctuations, at large A^c such 
disorder only occurs over distances exponential in A^c [26)127] . Neglecting such 
details, if chiral symmetry breaking occurs in vacuum, then it also occurs in 
the presence of a Fermi sea, in the following manner: 

($$)= cos(2/iz) ($'$') ; ($F5$) = -isin(2/i2) ($'$') . (39) 

This is the strict analogy of Migdal's pion condensation [22] in 1+1 dimensions. 
Chiral symmetry is broken, but by a spiral in the two possible directions, 
between $$ and $F^$. Schon and Thies [19] termed this as a "chiral spiral", 
and we adopt their evocative name, and so refer to our result as a Quarkyonic 
Chiral Spiral. 
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We can also understand why exciton pairing is not favored. An exciton conden- 
sate corresponds to which is a spin 1 operator. In the effective theory, 
the corresponding operators are $r^ri^2'^* and $r°ri 2$. These operators are 
flavor non-singlet. In two dimensions in vacuum, though, it is expected that 
the spontaneous breaking of chiral symmetry proceeds through condensates 
which are flavor singlets, and not through condensates which carry flavor. 

We conclude this section by noting that the extended flavor symmetry in 
1 + 1 dimensions is special to the vector-like interactions of QCD, and is not 
a generic property of theories in 3 + 1 dimensions. Consider, for example, 
a general Nambu-Jona-Lasino (NJL) model, with interactions such as {ipip)"^. 
This is the square of a mass term, and so is not invariant under the anomalous 



chiral transformation of Eq. (36 ). In accord with this. Nickel found that in NJL 
models at nonzero density, the thermodynamically favored ground state is a 
crystal, but not a chiral spiral [23]. 



5 Mapping Fluctuations for a Fermi Sea of Massive Quarks 



In this section we show how for massive quarks, excitations near the Fermi 
surface can be mapped onto a (modified) theory of the vacuum. For massless 
quarks, this could be done for the entire theory; here, it is only for fluctuations 
near the Fermi surface. 

Let the quark mass be m, so the Fermi momentum is related to the chemical 
potential as fJ''^ = p'p + rn?- We work in the extreme Quarkyonic limit, where 
/i ^ Aqcd- Starting with the theory in 3 + 1 dimensions, by neglecting the 
transverse momenta we obtain an effective theory in 1 + 1 dimensions. 



^T^'{^, + ^g2^A^)+^lT'^-nl^. (40) 



In this case, we could perform the transformation of Eq. (36), and so eliminate 
the term ~ yuF^ from the action, as in Eq. (37). The chemical potential is still 



in the action, though, through the transformation of the mass term into a 
complicated, position dependent "mass", as in Eq. (38). 



While one cannot make an exact correspondence to the vacuum theory, one 
can still make an interesting, if more limited, correspondence. Consider not 
all fluctuations in the theory, but just those near the Fermi surface. As is well 
known [Hj, even massive particles have a massless dispersion near the Fermi 
surface, with a modification to the speed of light: 



{pF + ^vY + m? — ^ Vf 5 p 



Vf 



\Pf 



(41) 
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As for the massless case, we can neglect anti-quarks, deep in the Fermi sea. 



First, let us decompose the spectrum of quarks. We have two Fermi surfaces, 
±|pp|, and have particle and hole excitations at each Fermi sea, with energy 
E ^ ±fir|(5p|. We identify the excitation around the Fermi sea with as a 

right moving fermion, and with — as a left moving fermion. The effective 
Lagrangian becomes 

^2d,massive ^^^^0 ^g^ ^ ^ ^ p. ^ ^ (42) 

where we have assumed = gauge. This is the same theory as for the 



massless case, Eq. (37), except the speed of light is not one, but vp. 



These elementary manipulations can be used to explain the nature of chiral 
spirals in exactly soluble models in 1 + 1 dimensions [T^PU] . 

For massless quarks, one maps the complete theory, in the presence of a Fermi 
sea, onto the vacuum. Thus the computation of ($'$') is complete, and both 
($$) and ($r^$) oscillate about zero. There are several examples, such as 
the Gross-Neveu model for a large number of flavors, in which this can be 
computed analytically [19]. 

Bringoltz considered a nonzero density of massive quarks for QCD in 1 + 1 
dimensions [20]. By numerical analysis of the theory in the canonical ensemble, 
he showed that chiral spirals also arise for massive quarks. These are due to 
quark and quark hole excitations about the edge of the Fermi surface, Eq. 



(42). As the density increases, vp — )■ 1, and the massive theory approaches the 



massless limit. 



6 Effective Theories for Excitations about the Fermi Surface 



The massless excitations near the Fermi surface can be described in terms of 
a Wess-Zumino-Novikov-Witten (WZNW) theory [2^26f27] . The dictionary 
between bosonic and fermionic currents for color, flavor, and (baryon) charge 
elements is 

=^tT[h-\d+h)tA] = , J- = ^tT[h{d-h-')tA] = :^ltA^-: , 

Ztt /vr 

tA is the color matrix, and h an element of SU{Nc); Tf is the flavor matrix, 
and g an element of SU{2Ni). Normal ordering of a composite operator A is 
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denoted by :A: 



For free Dirac fermions with colors and 2N{ flavors, 




(44) 



the bosonized version is 



S 



(45) 



with 



Sk[l] 




J (fx d^wn-^ + 



247r 



1 



/ 



£x e'"'\r^d^l){l-'dj){l-^dxl) . 



(46) 



This is a sum of a free massless scalar, for the U{1) of baryon number, a 
SU{Nc) WZNW model with level 2iVf, and a SU{2Ni) WZNW model with 
level iVc- 

A similar form can be derived for the theory of Dirac fermions coupled to 
a gauge field P^f^Tj . The flavor part of the action, 5^^™'' [(7], is completely 
unchanged, because the currents which define the flavor matrix g are color 
singlets. The color dependent part of the action, 5*2^°'' [/i], becomes that of a 
gauged WZNW action. 

The spectrum of a gauged WZNW model is involved [26|27] . However, what 
we are most concerned about are excitations near the Fermi surface; namely, 
are there gap less modes. In this context, what is of greatest concern is that is 
that the massless correlations of the flavor WZNW action, 5'^^™'' [(7] , dominate 
correlations over large distances. 

Note that this is true for both massless and massive quarks. In either case, 
there are numerous gapless modes about the Fermi surface. 



7 Conclusions 

In this paper we approximate the confining potential in a Quarkyonic phase 
as in Eq. We then find that Quarkyonic Chiral Spirals (QCS's) arise natu- 
rally, and can be expressed in terms of an effective model in 1 + 1 dimensions. 
Our analysis is valid in the extreme Quarkyonic limit, where /i 3> Aqcd, within 
a narrow skin of the surface of the Fermi sea, ~ Aqcd- 
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In this work we only considered the formation of a single chiral density wave, 
in a fixed direction. It is most likely, though, that the entire Fermi surface 
is covered with patches of chiral density waves, in different directions [TOIITT] . 
The detailed manner in which the Fermi surface is covered with such patches 
will be discussed separately [29] . 

Quarkyonic Chiral Spirals are reminiscent of the pion condensates of Migdal 
|22]- Pion condensates arise in effective models of nucleons interacting with 
pions: they are chiral spirals which oscillate as (■?/' exp(2i 0/^^75 where c is 
a pure number and /,r is the pion decay constant. In contrast, a QCS arises 
from the interactions of quarks and gluons; it is a chiral spiral not in chirality, 
75, but in spin, 7^7^. 

It is of great interest that we find a QCS for massive quarks. Kaplan and 
Nelson suggested that in dense matter, effective models of nucleons and kaons 
indicate that there is a condensate for the K~ field [23]. This is constant in 
space, and so is unlike the chiral spiral which we would expect in strange 
Quarkyonic matter. 

For distances x ^ QCS's have numerous modes with long ranged corre- 
lations, from the correlations of the WZNW model. Sec. |6| It is possible that 
the long range correlations of the modes of the WZNW model acquired finite 
range over distances ~ 1/Aqcd- Even so, this predicts man?/ more light modes 
than expected otherwise. 

The formation of QCS's has strong implications for the phase diagram of a 
gauge theory. The usual expectation is that one goes from a phase where 
chiral symmetry breaking is broken through a constant condensate, {ipip) 7^ 0, 
directly to a phase where it vanishes, {ipip) = 0. 

If a QCS forms, chiral symmetry is broken, but through an order parameter 
which differs from that in vacuum. As for a pion condensate, a QCS sponta- 
neously breaks both translational and rotational symmetries. Thus there is a 
strict order parameter which differentiates between ordinary chiral symmetry 
breaking, with a constant chiral order parameter, and a QCS. It also implies 
that there are exactly massless Goldstone bosons in a QCS, from the spon- 
taneously broken symmetries of translation and rotation. Similarly, there is 
an order parameter which differentiates between a QCS, and a phase which is 
chirally symmetric (at least for massless quarks). 

This implies that at zero temperature, as fi increases there is first a Fermi 
sea. There is then a well defined phase transition from a phase with constant 
{ipip) 7^ to a QCS. There is then a second phase transition, from the QCS to 
a chirally symmetry phase, with (ipilj) = 0. This differs from models which ex- 
clude spatially dependent condensates, Ref. [SIS], but similar to models which 
allow for then, Ref. [M] . 



21 



Thus at T = and 7^ 0, we predict the existence of an intermediate phase, 
with a QCS. In Appendix [A] we show that quarks affect gluons at asymptot- 
ically large /z ~ -^c^^- At such large fx, deconfinement and chiral symmetry 
breaking surely occur. Thus suggests that at least for large A^c, the region with 
a QCS is also large. Further, it is very possible that the phase transition from 
a QCS, to a chirally symmetry phase, occurs before the theory deconfines. 

We conclude by suggesting that dynamical quarks do not easily wash out 
a Quarkyonic phase, nor related effects, such as Quarkyonic chiral spirals. In 
QCD, it is known that the effects of screening, from dynamical quarks, are not 
strong. Notably, the linear term in the quark anti-quark potential persists to 
rather short distances, ~ 0.2 fm. This distance is comparable to the short range 
repulsion experienced by nucleons. Short distances then translates into high 
densities. It would be very interesting to estimate this effect within effective 
models, such as approximate solutions to the Schwinger- Dyson equations. 

Our analysis clearly raises more questions than it answers. However, we hope 
that we have provided a different way for thinking about cold, dense quark 
matter, and about the surprises which it might provide. In the end, this is of 
direct relevance to neutron (Quarkyonic?) stars, for which a wealth of exper- 
imental data will be forthcoming in the next few years. 
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A Screening by Dynamical Quarks at Large Nc 



In this appendix we sliow tliat at a nonzero temperature, T, screening by 
dynamical quarks enters when /i ~ N^^'^T. 

Consider, for example, the square of Debye screening mass at one loop order 
[6J. Gluons contribute ~ g'^N^.T'^, quarks give ~ g^Nijji^. Balancing the two 
terms, quarks start affecting gluons when when ~ N^^'^T. 

Next consider the same analysis for the free energy. Contributions from gluons 
are ~ T^, while those from quarks are ~ iV^A^f times powers of /i and T: /i'^. 



/i^T^, and T"^, see Eq. (A.3) below. Of course at ~ g'^ and beyond, the quark 
and gluon contributions get mixed up together, but this does not affect our 
power counting in A"c. If we then balance the leading term for gluons, ~ T^, 
against that of quarks, ~ NcN^fi^, we estimate that quarks dominate when 
/i ~ Nl^^T [6J, and not ~ N^^'^T, as for the Debye mass. 

With the free energy, however, one must take more care. One cannot simply 
equate the magnitude fo the free energies, but remember that deconfinement is 
only defined by a change in the relevant order parameter, which is the (renor- 
malized) Polyakov loop [30] ■ This is generated by a nontrivial distribution in 
the eigenvalues of the thermal Wilson line. This can be modeled by expand- 
ing about a constant expectation for the timelike component of the vector 
potential, 

Ao = ^. (A.l) 

9 

Q, is a diagonal matrix in color space, and is traceless, as a sum over elements 
of SU{Nc). Q has dimensions of mass, where the mass scale is set by the 
temperature, T. We do not need to know the explicit distribution of the Q's 
in order to estimate how large /i must be for quarks to affect the Q-distribution. 

At very large /x we can compute the quark determinant, in the presence of 
this background field, to one loop order. A constant field for Aq acts like an 
imaginary chemical potential for color, and so it is natural to introduce a color 
dependent chemical potential [3T] . 

Jl = fx + iQ. (A.2) 

At one loop order one obtains the usual result for the pressure, with replaced 
by yu: 

Pquark = A^f i — tr Jl^ + — ti Jl^ + — tM . (A.3) 

For large fi, the dominant contribution is from the expansion of the first term, 
~ Af tr /I*^ ~ fj,'^, as we estimated above. While this is as large as the gluon 



23 



contribution when /i ~ N^^^T, that does not matter, since this term is inde- 
pendent of Q. This holds order by order in perturbation theory, simply because 
Q has dimensions of mass. The next term is from the expansion of ~ A^f tr yU*^, 
equal to ~ n^tiQ] this vanishes, though, because Q is a traceless matrix. 

The leading term which is Q-dependent is the next term in the expansion of 
tr /i^, which is ~ /i^ tr Q^. Since the trace is ~ A^'c, this is as large as the gluon 
contribution, ~ N^, when /i ~ Ari/2T. This agrees with our estimate from the 
Debye mass. 

At large N^, since the quarks do not affect gluons until /i ~ N^/^T, the 
boundary from the confining, to the deconfining, phase is a straight line in 
the plane of T and fi [6]. The pure glue theory has a global symmetry of 
Z{Nc), which for large Nc is approximately U{1). At large /i, the leading term 
from quarks is ~ trQ^, and breaks the Z{Nc) symmetry. Like other terms 
from quarks, this favors a real expectation value for the Polyakov loop. Such 
a term acts to wash out the line of first order transitions. 

There are two possibilities. One is that the quarks produce a critical endpoint 
for deconfinement. The other is that the first line for deconfinement bends and 
meets the axis for T = 0. Since all of our analysis depends upon T 7^ 0, we 
favor the former. 



B Anomalous Baryon Number 



In this appendix we compute the baryon number generated by the anoma- 
lous transformation in Eq. (36). Consider the operator for baryon number, 
computed with point splitting: 



lini($(x + e)r°$(x))^^o = lim 



lim 

e-S>0 



=2 lim 



0^ _Z[fi ^ 0] 

Z[f^ = 0] 
tr[rV]+^/ie,tr[r^rV] + 0(e3) 



\ 27re2 



=2 hm 



(B.l) 



We take the symmetric limit to preserve Lorentz symmetry in the ultraviolet 
regime. 



lim 



(B.2) 
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The first term is odd in ~ eo, and so vanishes after averaging over both 
directions. Hence 

lini($(x + e)r°$(x))^^o = - . (B.3) 

This equals the baryon number density for a gas of two flavors of free quarks, 
with chemical potential /j,. 

The computation can be repeated for other operators, such as ($$), ($r^$), 
and ($r^$). These operators do not have anomalous terms, since 

tr[r5rV] = o, tr[rV] = o, Jtr[r^rv] = ^ ^ o , (b.4) 

respectively. For these condensates, the explicit phase dependence, in trans- 
forming from $ to must be taken into account, but there are no additional 
terms from the anomaly. 
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